Examining the Location of the Magnetopause in an Undergraduate Lab 

James P. Crumley, Ari D. Palczewski, and Stephen A. Kaster 
College of Saint Benedict / St. John 's University Collegeville, Minnesota Q 

(Dated: June 13, 2007) 

Abstract 

Integrating areas of current research into undergraduate physics labs can be a difficult task. The location 
of the magnetopause is one problem that can be examined with no prior exposure to space physics. The 
magnetopause location can be viewed as a pressure balance between the dynamic pressure of the solar wind 
and the magnetic pressure of the magnetosphere. In this lab students examine the magnetopause location 
using simulation results from BAT-R-US global MHD code run at NASA's Community Coordinated Model- 
ing Center Students also analyze data from several spacecraft to find magnetopause crossings. The students 
get reasonable results from this lab as well as exposure to the tools and techniques of space physics. 
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I. INTRODUCTION 



Physics laboratory courses often emphasize experiments that relate directly to physical con- 
cepts being covered in lecture courses that students are taking concurrently. These experiments 
are intended to allow students to see for themselves how physical concepts work in the real world. 
Experiments of this sort tend to deal with physics that has been understood for some time and of- 
ten involve performing classic experiments of the past. Experiments of this type do add to student 
understanding, but a steady diet of them can leave students disconnected from the current practice 
of physics. Maintaining the students' intellectual curiosity about science is a key to retaining stu- 
dents in science majors.^ One way to maintain that curiosity is to create experiments that expose 
students to areas of current research to go along with classical experiments. 

Designing experiments that deal with areas of current research and are accessible to undergrad- 
uate students can be a difficult task. This task is particularly difficult for space physics because 
most students have had little exposure to space and plasma physics. Space physics applications 
are often complex, defying simple treatment, and they often rely on advanced electricity and mag- 
netism, which students often have late in their undergraduate curriculum. 

In this paper we will discuss our use of the location of Earth's magnetopause as a topic for a lab 
for physics majors. The magnetopause is defined as the boundary between Earth's magnetic field 
and the interplanetary magnetic fields. The magnetopause is an appealing topic for a lab because it 
can be introduced at a fairly elementary level as a pressure balance, after which more complicated 
models can also be examined. This lab involves examining magnetopause location using computer 
simulations and spacecraft data. 

In Sec. In] we introduce some background information regarding the magnetosphere in general, 
and the magnetopause in particular. The portion of the lab using computer simulation to examine 
the magnetopause location is discussed in Sec. Hill while the portion of the lab using spacecraft 
observations is discussed in Sec. |IVl Conclusions as well as ideas about similar sorts of labs are 
examined in Sec.lVl 
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II. MAGNETOSPHERE INTRODUCTION 



A. The Magnetosphere and the Magnetopause 

The structure and behavior of the area of space just outside of the Earth's neutral atmosphere is 
a problem that physicists have worked on for quite some time. The first clues to the nature of this 
region, now known as the magnetosphere, came from observations of the magnetic field on the 
Earth and observations of comets.^ Though the Earth's magnetic field had been used since ancient 
times for navigation, in 1600 William Gilbert was the first to purpose that the Earth was a giant 
magnet.^ Serious observations of the Earth's magnetic field in the 1700s led to the discovery of 
variations in the magnetic field of the Earth called magnetic storms."* Magnetic storms were long 
hypothesized to be due to the Sun.- It took some time to arrive at a suitable mechanism^ for the 
Sun to be affecting the magnetosphere, but it is now known that the solar wind causes magnetic 
storms. Comets were seen to have two tails: one caused by light from the Sun, the other by a 
stream of particles from the Sun, now called the solar wind.'^ 

The magnetic field of the Earth acts as an obstruction to the solar wind, causing a shock called 
the bow shock and a boundary region called the magnetosheath (see Figured]). The interaction 
with the solar wind compresses the Earth's magnetic field on the side of the Earth toward the 
Sun and stretches it on the side away from the Sun, resulting in the asymmetrical shape of the 
magnetosphere. As a first approximation, the location of the magnetopause is set by the balance 
between the dynamic pressure of the solar wind and the magnetic pressure of the magnetosphere. 
Variations in the solar wind lead to movements in where the magnetopause is located, as well 
as variations in the magnetic field measured on Earth. Currents within the magnetosphere and 
ionosphere also affect the magnetic field measured on Earth.- 

The existence of a magnetosphere with a magnetopause of this sort was first posited to help 
explain variations in the readings of Earth-bond magnetometers^ and to tie those variations to 
processes on the Sun. After these predictions, early spacecraft missions confirmed the existence 
of both the magnetosphere and the magnetopause.'' Since that time there has been advancement in 
the observation, theory, and modeling of the location of the magnetopause. Recent models attempt 
to fit spacecraft observations of the magnetopause location and solar wind conditions to a variety 
of functional forms. 
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B. Magnetopause Pressure Balance 



As mentioned above, the location of the magnetopause can be considered to be due to the 
pressure balance between the dynamic pressure of the solar wind and the magnetic pressure of 
Earth (see Figure |2l): 

'2psw{Vsw COS ef = ^^B^nside- (1) 

where psw is the mass density of the solar wind, Binside is the magnetic field inside the magne- 
topause, and 6 is the angle of incidence of the solar wind (the v^^ cos 6 term is the solar wind speed 
normal to the magnetopause). 

It turns out that due to the compression of the magnetosphere Binside is not simply the dipole 
magnetic field of the Earth: 

B dipole = Boi^^—^-^ , (2) 

where Bo is the surface magnetic field at the equator, i?^; is the radius of the Earth, and r is 
the distance from the center of the Earth to the location of interest. In this Chapman-Ferraro 
Model'° for the magnetopause boundary location, there is a current that runs tangential to the 
magnetopause boundary in the dawn to dusk direction. This current causes a magnetic field that 
cancels the magnetic field of the Earth just outside the magnetopause, and doubles the magnetic 
field inside the magnetopause^ so that: 

Binside '^Bdipole 2i?Q I 1 . (3) 

\ r J 

Combining Eqs.[il[2l and[3]and solving for the location of the magnetopause leads to: 

r ( Bl \^ 



(4) 



Re KPoPswivswCos^ 6) 
Assuming that solar wind consists of protons and electrons coming in normal to the magnetosphere 
and substituting in typical values of 10 protons/cm^ and 400 km/s, Eq. |4] gives a distance to the 
magnetopause subsolar point, r^, of roughly 8 Re- (The subsolar point is the location on the 
magnetopause along the line from the Earth to the Sun.) The observed value is for those conditions 
is 10. Re. 

A more thorough empirical expression including several factors ignored above and taking ^ = 

is 

r„(i?E) = 107.4 (n,^t;L)"^. (5) 
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In Eq. [51 To is the distance from the center of Earth to the magnetopause subpolar point in Rg, 
Usw is the number density of the plasma in the solar wind in cm^'^, and Vsw is the speed of the 
solar wind in km/s.— Note that this expression has the same dependence on solar wind speed and 
number density as Eq. H but a different leading constant. 

C. Recent Magnetopause Modeling 

Much work has been done on modeling the location and shape of the magnetopause since 
Chapman and Ferraro's original model.— Along with the dynamic pressure of the solar wind, it 
was found that the orientation of the interplanetary magnetic field plays a key role in determining 
the shape of the magnetopause.— More recent studies have concentrated on using databases of 
magnetopause crossings of various spacecraft and the solar wind conditions at those times to 
formulate empirical expressions for the magnetopause location.^^^ In the lab described here we 
use fits from the work of Shue et al.-^^^ 



In those Eqs. [6]-[8l r is the distance from Earth to the magnetopause boundary in Re, is the 
distance from Earth to the subsolar point of the magnetopause in Re, is the z-component of 
the solar wind's magnetic field in nT, Dp is the dynamic pressure of the solar wind in nPa, and a 
is a unitless number representing the amount of tail flaring on the night side of the magnetosphere. 

ni. MAGNETOPAUSE LOCATION FROM SIMULATION 

In this portion of the lab students perform computer simulations using various solar wind con- 
ditions and determine where the subsolar point of the magnetopause is from their results. These 
results are fit to an expression of the form of Eq. [51 with the leading constant as a free parameter. 
The students compare their constant to the value of 107.4 from Eq.[51 




(6) 



(10.22 + 1.29 tanh [0.184(5^ + 8.14)])(Dp)^ 
a = (0.58 - 0.007S^) [1 + 0.024 In (Dp)] 



(7) 



(8) 
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A. Simulation Environment 

The students model the magnetosphere by running the BAT-R-US^ simulation on supercom- 
puters at the Community Coordinate Modeling Center dCCMCl) . BAT-R-US solves the 3D mag- 
netohydrodynamics equations in finite volume form using an adaptive grid. The simulation is run 
on NASA supercomputers at the CCMC where simulation runs can be requested by the public. 
Results from the simulations can be explored and visualized through the use of a standard web 
browser. The documentation provided is quite thorough, so explaining to students how to request 
a run and access the results is relatively easy. The resources that CCMC provides allow for access 
to powerful simulations without the need to have local access to supercomputer hardware or to 
install and maintain simulation codes. 

B. Finding tlie Magnetopause Location in Simulation Results 

In the process of finding the subsolar point of the magnetopause in their simulation results, each 
lab group was required to come up with their own standards for determining the magnetopause lo- 
cation. The lab materials explain what processes will be going on near the magnetopause and how 
that might affect the students' graphs, but no prescription for finding the magnetopause is given. 
The task of defining their own standards serves several purposes. The students are stimulated by 
this requirement to examine their graphs more closely than they would be if a prescriptive method 
for finding the magnetopause were given. Furthermore, the openness of this task is a fair reflection 
of real lab work since interpretation of simulation results often requires discretion. 

The instructions regarding finding the magnetopause location advise the student to look at more 
than one plasma parameter when finding the magnetopause, since the magnetopause should be 
evident in more than one type of data. Since the magnetopause is defined as the boundary between 
where the sun and the earth's magnetic fields dominate, it is clear that a signal should be visible 
in the magnetic field results. The magnetopause also affects the motion of the plasma particles, so 
differences in the plasma bulk velocity and number density should also be visible. 

Figure [3] shows an example of data used by students to find the subsolar point of the magne- 
topause. On this plot, the magnetopause is located at roughly 11 R^. In the number density plot, 
the bump in the number density corresponds to the buildup of plasma in the magnetosheath, so the 
inner boundary of that bump corresponds to the magnetopause. In the plot of the x-component of 
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the velocity, the magnetopause is seen as the location where the value goes to 0, since the plasma 
from the solar wind is diverted around the magnetosphere at the magnetopause. Finally, in the plot 
of the z-component of the magnetic field, there is an almost imperceptible shift at 11 R^;. So in 
this case, counterintuitively, it is easier to find the magnetopause in the plasma results for the the 
simulation than in the magnetic field results. 

C. Empirical Fit 

In this lab, students find the subsolar point of the magnetopause for a variety of plasma con- 
ditions and fit their data to Eq. [51 finding their own value to compare to the leading constant of 
107.4. The students run the simulations for two types of conditions. In one set, the solar wind 
speed is held constant and the solar wind number density is varied, and in the other set the fixed 
and varying parameters are switched. Most groups get good agreement between their leading con- 
stants for these two sub-groups. Some groups have had constants that were a bit high. Finding 
high values for the constant is not surprising since the students were given so much freedom to 
define their own method of determining the magnetopause location. More prescriptive directions 
on how to find the magnetopause location might yield better constants, but we are wary of losing 
the experience that the students gain by defining their own standards. 

rV. SPACECRAFT DATA 

In this section of the lab, students search for magnetopause crossings for three sets of data cho- 
sen from several events and several spacecraft (Geotail, Polar, the GEOS satellites, and the LANL 
geosynchronous satellites). The students first compare the actual position of the spacecraft to the 
location of the magnetopause predicted by Eq. [6l Next, they examine the particle and magnetic 
field data from that spacecraft for signs that the spacecraft made magnetopause crossings at the 
predicted time or at other times. Finally, students compare the signs of magnetopause crossings 
that they see in the magnetic field and particle and discuss their perceptions of the difficulty of 
interpreting those two types of data. 
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A. Magnetopause Crossing in Data 



Figure m shows an example of the search for spacecraft crossings of the magnetopause, in this 
case Geotail, that students do in this lab. The data shown is from 31 October 2003. On that day a 
large CME hit Earth, causing auroras that were visible throughout much of the United States.— 

The top panel of Figure |4] predicts several magnetopause crossings, but the most notable cross- 
ings are at roughly 5:00, 10:00, and 11:00. From 1:30-5:00 and from roughly 10:00-11:00, Geo- 
tail is predicted to be inside the magneto sphere. In other words, at those times Figure |4] shows the 
magnetopause further from the Earth than the spacecraft is. For most of the rest of that day Geotail 
was inside the magneto sphere. The data in the middle and bottom panels of Figure |4] shows the 
crossing from inside to outside the magnetosphere at 1 1:00 and the data also shows other magne- 
topause crossings, though as described below the crossings seen in the data differ a bit from the 
predictions. 

The middle panel of Figure |4l which shows the ion flow velocity measured by Geotail and 
the x-component of the solar wind speed measured by ACE, has three visible regimes. During 
this day the x-component of the solar wind velocity slowly changes from —1200 km/s to —800 
km/s. From 0:00-5:00 most of the Geotail ion flow velocity data is missing. From 5:00-11:00 
the three components of the ion flow velocity measured by Geotail oscillate near values of —700 
km/s, —400 km/s, and 200 km/s respectively. While from 11:00-24:00 the y- and z-components 
oscillate near km/s, and the x-component is based at roughly —1000 km/s, with spikes up to 
almost km/s. The behavior of the x-component during this last time period can be interpreted 
as being due to the spacecraft being just outside the magnetopause in the magneto sheath. When 
the x-component is near —1000 km/s it matches the solar wind speed, suggesting that Geotail is 
outside the magnetopause at those times. The noisiness of these measurements of the ion speed in 
the magnetosheath is likely due to reflection of the solar wind ions off the magnetopause and the 
movement of the magnetopause as the solar wind conditions vary. 

The velocity measurements act much differently from 5:00-11:00, suggesting that the space- 
craft is inside the magnetosphere during those times. Notice that during this time there are several 
short time periods where all components of the data have spike which match their magnetosheath 
values. This suggests that during this time period the spacecraft is near the magnetopause bound- 
ary and that fluctuations in the solar wind cause the magnetopause to oscillate back and forth 
across Geotail's position. 
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The bottom panel, showing Geotail's magnetic field measurements, has similar regions of be- 
havior. From 0:00-1 1:00 all three components have spiky measurements which trend downward. 
The downward trend is due to the decreasing magnitude of the Earth's dipole field as Geotail gets 
further from Earth, while the spikes are due to disturbances in the plasma in the magnetosphere 
during this magnetic storm. From 11:00-24:00, the three components are steadier and oscillate 
near 10 nT, 5 nT, and —10 nT, respectively. During this time Geotail is in the magneto sheath where 
the solar wind magnetic field dominates. Note that during several of the spikes in the magnetic 
field which occur between 0:00 and 11:00, that the magnetic field matches the magnetosheath val- 
ues. This supports the idea that the magnetopause oscillated back and forth past the spacecraft, as 
was mentioned with the ion velocity data. 

In summary, the spacecraft data (bottom two panels of Figure U) show that Geotail crossed 
the magnetopause. The data suggest that Geotail was inside the magnetopause from 0:00-11:00 
(though there were some brief crossings during that time), and that it was outside the magnetopause 
from 1 1 :00-24:00. This is a slightly different story that what the top panel of of Figure|4]) predicted. 
The prediction had a more complicated series of crossings up until a final crossing to the outside 
of the magnetopause at 1 1:00. Overall though, the predicted and actual crossings agree fairly well. 

B. Results 

The students get reasonable results for this section of the lab. As seen in the example in 
Sec. IIV Al finding the magnetopause crossings in the data can be complicated process, so the 
students do typical need a bit of guidance as they proceed. One of the things that troubles some of 
the students is that some of the data sets they examine contain no magnetopause crossings. This 
disconcerts students since their expectation is that all data sets will have crossings. In research 
there are often data sets that do not contain the phenomenon being searched for, so it is good to 
have students do some cases of this sort. 

In general, the students do a good job of finding the crossings in the magnetic field data, but 
they have more difficulty with the particle data. These problems are to be expected since in most 
of the cases the ion data is not a clear as in Figure HI This problem could be alleviated somewhat 
by finding events where the signs of the magnetopause crossing are clearer in the ion data, as well 
as by giving the students more guidance in the interpretation of the ion data. 
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V. CONCLUSION 



We have presented a lab in which students examine the Earth's magnetopause using simula- 
tions and spacecraft data. In this lab students are challenged by exposure to an active area of 
research that they might not otherwise encounter as undergraduates. Furthermore, not only are 
they exposed to space physics, they work on a problem using real data and tools in a manner that 
is not far removed from what current researchers do. Experiences of this sort are key in retaining 
physics majors and in helping physics students determine what they want to do with their physics 
education. 

The outlook for future space physics and astronomy labs involving nearly current research 
topics is bright. The move toward open access to spacecraft data, coupled with the increase in the 
number of operational spacecraft taking science data, continues to broaden the areas of space that 
anyone with an internet connection can have access to. On the simulation side, the development of 
efforts to give public access to supercomputers, along with the steady increase in computing power, 
continues to expand the number of problems that anyone can find the resources to simulate. The 
key difficulty to applying these new observational and computational powers to the teaching lab 
remains formulating appropriate problems for students to examine. 
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Figure Captions 




FIG. 1 : The Earth's magnetosphere. 




Magnetopause 



FIG. 2: The pressure balance at the magnetopause between the solar wind dynamic pressure and the Earth's 
magnetic pressure. 
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FIG. 3: Simulation results plot used to find subsolar point of the magnetopause. The data are from five 
minutes into a typical simulation run and they are plotted along the line from the Earth to the Sun. The 
top panel shows the number density, the middle panel shows the x-component of the plasma flow velocity, 
and the bottom panel shows z-component of the magnetic field. Geocentric solar magnetospheric (GSM) 
coordinates are used here. In this system the x-axis points from the Earth to the Sun, the z-axis points in the 
direction of Earth's north magnetic pole, and the y-axis completes the right-handed system.— 
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FIG. 4: Plots showing Geotail's magnetopause crossings on 31 October 2003 (color online). The top panel 
shows the Geotail spacecraft's distance from Earth and the predicted location of the magnetopause along 
the line from the Earth to Geotail based on the ACE spacecraft measurements of solar wind conditions. 
The middle panel shows Geotail's measurements of the ion flow velocity and the x-component of the solar 
wind speed measured by ACE, all in geocentric solar ecliptic (GSE) coordinates. In this system the x-axis 
points from the Earth to the Sun, the z-axis points perpendicular to the plane of earth's orbit, and the y-axis 
completes the right-handed system.^** The bottom panel shows the three components of the magnetic field 
measured by Geotail in GSE coordinates. Note that 100 nT has been adding to By in order separate the hnes 
on the plot. 
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